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' the evolution of primordial density perturbations. If the phase transi- 

ts . tion is first order, the sound speed vanishes during the transition and 

density perturbations fall freely. The Harrison-Zel'dovich spectrum 
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Abstract 

We investigate the effect of the cosmological QCD transition on 



of density fluctuations develops large peaks on scales below the Hub- 
ble radius at the transition. These peaks above the primordial spec- 
trum grow with wavenumber and produce cold dark matter clumps of 
£2 \ masses less than 1O -1O M0. At the horizon scale the amplification of 

overdensities is of order unity, and thus no 1M black hole production 
is possible for a COBE normalized scale-invariant spectrum. 



Introduction. The QCD transition from a hot plasma of deconfined quarks 
and gluons to a hot gas of hadrons happens in the early Universe at a tem- 
perature T*. Lattice QCD results for the physical values of the quark masses 
indicate that the QCD phase transition is of first order |] and takes place 
at T* ~ 150 MeV ||. The QCD transition may have important cosmological 
consequences: it could lead to inhomogeneous nucleosynthesis ||, produce 
dark matter in the form of strangelets |3| , or give rise to gravitational waves. 
These effects have in common that their scale is set by the mean bubble 
nucleation distance -R nU ci and not by the Hubble scale i?n- It has turned out 
that -R nU ci <C -Rh, mainly because the surface tension is very small ||. In 
contrast, we study the scales A 3> -R nU ci- 

Cosmological density perturbations are affected by the QCD transition 
at scales A ~ i?n H 0, where i?n ~ m p/T* ~ 10 4 m. For a first order 

*Talk given at the Third International Conference on Physics and Astrophysics of the 
Quark Gluon Plasma, March 17 - 21, 1997, Jaipur, India. 
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QCD phase transition the deconfined and confined phases can coexist at 
the coexistence temperature T* at fixed pressure = p(T+). During this 
coexistence period pressure gradients, and thus the sound speed, vanish for 
wavelengths A ^> R nuc i- The vanishing sound speed gives rise to large peaks 
and dips above the primordial spectrum of density perturbations. These 
peaks grow at most linearly with wavenumber. We show that the formation 
of black holes at the horizon scale is impossible for standard inflationary 
scenarios, in contrast to recent claims [^, pi- We predict clumps with M < 
l(r lo M in cold dark matter (CDM), if the CDM is kinetically decoupled at 
the QCD transition. 

The cosmological QCD transition. The QCD phase transition starts 
with a short period (10~ 4 £h) of tiny supercooling, 1 — T sc /T* ~ 10 -3 . When 
T reaches T sc , bubbles nucleate at mean distances -R nuc i ~ 2 cm |]TT|. The 



bubbles grow most probably by weak deflagration ||10| . The released energy 
is transported into the deconfined phase by shock waves, which reheat the 
deconfined phase to T* within 10~ 6 tH- The entropy production during this 
very short period is dS/S ~ 10^ 3 . Further bubble nucleation in the reheated 
quark-gluon phase is prohibited after this first 10 _4 tn- Thereafter, the bub- 
bles grow adiabatically due to the expansion of the Universe. The transition 
completes after lO -1 ^- 

During this equilibrium evolution the rate at which the quark-gluon phase 
can be converted into the hadron phase and vice versa, r conv , must be com- 
pared to the Hubble rate H and to the wave number fc phys of the considered 
density perturbation. r conv cannot enormously differ from the typical QCD 
scale fm -1 , while H « (10 4 m) _1 , therefore H/T conv and fc phys /r conv are of 
order 10~ 19 . This means that the QCD relaxation times are incredibly short 
and that the 2-phase system is very close to equilibrium. Since all interaction 
rates r strongi eiweak H, all particles are in chemical and thermal equilibrium. 
At scales A > -R nuc i we may describe photons, leptons, and the QCD matter 
by a single radiation fluid with equation of state p = p(p). 

During this reversible coexistence period (isentropic) T and p are fixed, 
while p can vary, therefore 

/ / in \ 

= 0. (1) 

isentropic 




In [0 it was claimed that the isentropic condition is not consistent. This is 
wrong, because they apply non-relativistic (u < l,p < p) hydrodynamical 
equations to relativistic wave fronts in the radiation fluid. Moreover, they 
claim that sound propagation keeps the energy fraction in each phase fixed. 
This is incompatible with the equilibrium situation described above. 
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Figure 1: The modifications of the CDM density contrast |<5cdm|(^*/10) and 
of the radiation fluid amplitude Arad due to the QCD transition (lattice 
QCD fit). Both quantities are normalized to the pure Harrison-Zel'dovich 
radiation amplitude. On the horizontal axis the wavenumber k is represented 
by the CDM mass contained in a sphere of radius n/k. 

The thermodynamics of the cosmological QCD transition can be studied 
in lattice QCD, because the baryochemical potential is negligible in the early 
Universe, i.e. ph/T* ~ 10~ 8 . We fit the equation of state from lattice QCD 
|§. In addition we use the bag model to illustrate the origin of the large 
peaks in the density spectrum. 

The evolution of density perturbations. During an inflationary epoch in 
the early Universe density perturbations 5p have been generated with a scale- 
invariant Harrison-Zel'dovich (HZ) spectrum. During the radiation domi- 
nated era subhorizon (k phys 3> H) density perturbations oscillate as acoustic 
waves, supported by the pressure 5p = c 2 5p. It is useful to study the di- 
mensionless density contrast 5 = 5p/p as a function of conformal time and 
comoving wave number. The subhorizon equation of motion reads 

S" + c 2 s k 2 5 = , (2) 

if we assume that the phase transition is short compared to the Hubble time. 
From the lattice QCD equation of state the period of vanishing c 2 s lasts for 
O.ltn, in the bag model it lasts for 0. Bin- 
To illustrate the mechanism of generating large peaks above the HZ spec- 
trum we discuss the subhorizon evolution of the density contrast S in the bag 
model J7[. Before and after the QCD transition c 2 = 1/3 and the density con- 
trast oscillates with constant amplitudes A m resp. A out . At the coexistence 
temperature T* the pressure gradient, i.e. the restoring force in in Eq. (^), 
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vanishes and the sound speed ([!]) drops to zero. Thus 5 will grow or decrease 
linearly, depending on the fluid velocity at the moment when the phase tran- 
sition starts. This produces large peaks in the spectrum which grow linearly 
with the wave number, i.e. A ont / A m \ pea x iS ~ k/kx- The radiation in a Hubble 
volume at T* has the mass Mrad(-Rh) ~ 1M . The CDM mass is related to 
the radiation mass by Mcdm = a,*/a eq MRAD ~ 10 -8 Mrad- CDM falls into 
the gravitational potential wells of the radiation fluid during the coexistence 
regime. 

Fig. H shows the final spectrum for the density contrast of the radiation 
fluid and of CDM. It has been obtained from numerical integration of the 
fully general relativistic equations of motion and the lattice QCD equation of 
state . The mass Mi corresponds to the wave number k\ of the bag model 
analysis. M\ coincides with the mass inside the Hubble horizon. The mass 
M 2 stems from a WKB analysis with the lattice QCD equation of state, which 
shows that A ont / A m \ peaks « (k/k 2 ) 3 ^- We conclude that at the horizon scale 
the modification of the HZ spectrum is mild, whereas for scales much smaller 
than the horizon big amplifications are predicted. Without tilt in the COBE 
normalized spectrum the density contrast grows nonlinear for k phys /H ~ 10 4 
resp. 10 6 for the bag model resp. lattice QCD equation of state. 

Implications of the large peaks. It was suggested that the QCD transi- 
tion could lead to the formation of 1M black holes that could account for 
all missing (dark) matter || [|, i.e. all the radiation mass within one Hubble 
horizon should collapse to a single black hole. In common inflationary sce- 
narios the density contrast is of order 10~ 4 on subhorizon scales before the 
QCD transition. Thus, our linear analysis applies. From Fig. p] it is clear 
that there is only a small amplification of order unity at the horizon scale 
(around Mi). Therefore, with COBE normalized scale-invariant spectrum 
(allowing tilts \n — 1| < 0.3) it is impossible to form black holes at the QCD 
transition that are abundant enough to close the Universe. 

At T ~ 1 MeV the neutrinos decouple from the Hubble scale, from smaller 
scales they decouple at higher temperatures. During their decoupling they 
damp all fluctuations in the density contrast due to collisional damping. 
Thus, at the time of big bang nucleosynthesis all peaks of Fig. |] are erased 
and the energy density on scales A ~ Rn(l MeV) is homogeneous. 

Weakly interacting dark matter, like the neutralino, does not belong to 
our CDM, because it is not kinetically decoupled at the QCD transition. For 
our purpose the neutralino would belong to the radiation fluid and subhorizon 
perturbations in its density are damped during its kinetic decoupling. 

However, collisional damping is irrelevant for other CDM like primordial 
black holes. Peaks in this CDM survive and grow logarithmically during the 
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radiation era. Shortly after equality they grow nonlinear and collapse by 
gravitational virialization to clumps of M < 10~ 10 M Q . Their size can be 
estimated to be a few AU. Whether such CDM clumps may be observable 
remains to be investigated. 

D. J. S. would like to thank the organizers of this conference for making 
stimulating talks and discussions possible. D. J. S. and P. W. are supported 
by the Swiss National Science Foundation. 
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